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Abstract 

We consider the observational constraints from the detection of antiprotons 
in the Galaxy on the amount of Primordial Black Holes (PBH) produced 
from primordial power spectra with a bumpy mass variance. Though es- 
sentially equivalent at the present time to the constraints from the diffuse 
7-ray background, they allow a widely independent approach and they 
should improve sensibly in the nearby future. We discuss the resulting con- 
straints on inflationary parameters using a Broken Scale Invariance (BSI) 
model as a concrete example. 

PACS Numbers: 04.62. +v, 98.80. Cq 



1 Introduction 



The formation of PBHs in the early universe is an inevitable prediction based 
on general relativity, the existence of a hot phase and, most importantly, the 
presence of primordial fluctuations which are the seed of the large structures in 
our universe It can have many interesting cosmological consequences and is 
one of the few constraints available on the primordial fluctuations on very small 
scales that can be based on existent astrophysical observations (see e.g. 0). It 
has been used by various authors in order to constrain the spectrum of primordial 
fluctuations, in particular in order to find an upper limit on the spectral index 
n and on the present relative density of PBHs with M (the initial mass 

of a PBH whose lifetime equals the age of the Universe) A possible 

contribution of evaporating PBHs to the diffuse 7-ray background is presently 
the most constraining observation 0. On the other hand, the observation of 
antiprotons in the Galaxy is as powerful § and, in contrast to the 7-ray 
background, sensitive improvements can be expected in the near future. These 
involve both experimental and theoretical progress. This is why it is interesting 
to consider in some details the constraints these observations can, and will, put on 
any primordial fluctuations model, and prominently on some inflationary models. 
As noted earlier (see, e.g., Fig.l in 0), a constant spectral index n would need 
extreme fine tuning in order to saturate the 7-ray or antiproton constraint, and 
such a large n is anyway excluded by the latest CMB data. Hence we consider here 
spectra with a characteristic scale for which the generation of PBHs is boosted 
in a certain mass range. 



2 PBH formation and primordial fluctuations 

Density of PBHs from humpy mass variance: For detailed confrontation with 
cosmological and astrophysical observations one often needs the mass spectrum, 
the number density per unit of mass. This is particularly delicate for PBHs and 



we follow here a derivation valid in the presence of a bump, as given in |]T0 
The first assumption is that the primordial spectrum of cosmological fluctuations 
has a characteristic scale in its power spectrum P{k), which results in a well- 
localized bump in its mass variance. The importance of this assumption lies 
in the determination of the PBH mass scale Mp^ak where PBH formation mainly 
occurs. The second assumption, supported by numerical simulations, is that PBH 
formation occurs through near-critical collapse |[TT[| whereby PBH with different 
masses M around Mpeak = Mnitk^,^^,,), the horizon mass at the (horizon-crossing) 
time tfcpeafc ^ horizon crossing time is defined through k = a(tk)H{tk) - 
could be formed at the same time tu , , according to 

"'peak ' ^ 

M = kMh{6-6,)\ (1) 

where 6c is a control parameter. While the parameters 7 and 6c are universal 
with 7 !^ 0.35, 6c ~ 0.7, the parameter k, (or e, see below) can vary sensibly and 
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fixes essentially the typical PBH mass. As shown in |]TU[, one finds 



dQpBH _dnpBHiM,t,J / / M 
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If we identify the maximum of (j^) in the following way 

Mmax = e Mpeak , (3) 

we are led to the result 

-e~^(l + 7) 



e ^ /3{Mpeak) (1 + -) TT exp 



'-peak 



M \ ' 



peak 



(3{Mpeak) gives the probability that a region of comoving size R = {H ^ j a)\t=t 
has an averaged density contrast at the time tfop^a* range 5c < < ^max 

^max 



(4) 



peak 



mveak)= piS^h^^JdS. (5) 

It is then straightforward to find the quantity of interest to us 
(Pm / Mp \^ _2 dVtpBH 



X — — —{x) , (6) 



dMidVi 327r \ Mpeak J dlnM 

where M„ stands for the Planck mass while x = .J^ . The subscript i stands 
for "initial", i.e. at the time of formation. The mass Mp^ak corresponds to the 
maximum in the mass variance anitk) and not to the maximum in the primordial 
spectrum itself 0. The parameters 7 and e refer to PBH formation while Mpeak 
and P {Mpeak) refer to the primordial spectrum. 



Primordial inflationary fluctuations: One usually considers Gaussian primordial 
infiationary fiuctuations but it should be stressed that non-Gaussianity of the 



fiuctuations could lead to sensibly different results |T2[- For primordial fluctua- 



tions with a Gaussian probability density p[6], we have 

1 1 r"^ 

p{5) = , a^R) = — dk W^^{kR) P{k) , (7) 

V27r a(R) 2tt^ Jq 



27r a{R) 

where 5 is the density contrast averaged over a sphere of radius i?, and = 
^ j ^ is computed using a top-hat window function. Usually what is meant 

by the primordial power spectrum is the power spectrum on superhorizon scales 
after the end of inflation. On these scales, the scale dependence of the power 
spectrum is unaffected by cosmic evolution. On subhorizon scales, however, this 
is no longer the case, and one has instead 

P{K t) = Pik, t,) T'ik, t) , T(fc 0, t) 1 , (8) 
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where ti is some initial time when all scales are outside the Hubble radius {k < 
aH). Therefore, the power spectrum P{k) on sub- horizon scales appearing in 
must involve convolution with the transfer function at time tjt 0]. At reen- 
trance inside the Hubble radius during the radiation dominated stage, one has in 
complete generality [13], (subscript e stands for the end of inflation) 



8l7r2 



F{kx) T^{kx,tk) W^Hi^) da; 



tfc, < tfc < teg , (9) 



where the transfer function can be computed analytically and yields 



T'^{kx,tk) 



9 / sin(csx) 
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while F{k) = %k^P{k,tk) 



-iT'^6'fj{k,tk). Finally (3{Mpeak) is given by 



(Tnit) 



peak ) 
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and we will take Sr = 0.7. 



with aUtk^^J ^ ^\R)\t 

For a given primordial fluctuations spectrum of inflationary origin normal- 
ized at large scales using the COBE data, the quantities Mp^ak and (3{Mpeak) 
can be computed numerically and will depend on some inflationary parameters 
specifying that model as well as on cosmological parameters pertaining to the 
cosmological background evolution [l^. On the other hand 7 and e should be 
found by numerical simulations of PBH formation for this particular spectrum. 
Values e = 0.5, 1, 2, correspond to k ^ 2.7, 5.4, 10.8. 



3 Evaporation, fragmentation and source term 

As shown by Hawking , such PBHs should evaporate into particles of energy 
Q per unit of time t (for each degree of freedom): 

(12) 



h (exp (^) - (-1)- 

where contributions of angular velocity and electric potential have been neglected 
since the black hole discharges and finishes its rotation much faster than it evap- 



orates [|T5|. The quantity k is the surface gravity, s is the spin of the emitted 
species and F^ is the absorption probability. If the Hawking temperature, defined 
by T = hc^ / {IQukG M) ^ (lO^^g/M) GeV is introduced, the argument of the ex- 
ponent becomes simply a function of Q/kT. Although the absorption probability 
is often approximated by its relativistic limit, we took into account in this work 
its real expression for non-relativistic particles: 
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where as is the absorption cross section computed numerically |T^ and /i is the 
rest mass of the emitted particle. 

Among other cosmic rays emitted by evaporating PBHs, antiprotons are es- 
pecially interesting as their secondary flux is both rather small (the p/p ratio 
near the Earth is lower than 10~^ at all energies) and quite well known [Q. We 
will therefore focus on such antiparticles in this paper. As shown by MacGibbon 
and Webber when the black hole temperature is greater than the quantum 



chromodynamics confinement scale Aqqd, quarks and gluons jets are emitted in- 
stead of composite hadrons. To evaluate the number of emitted antiprotons p , 
one therefore needs to perform the following convolution: 



where aj is the number of degrees of freedom, E is the antiproton energy and 
dgjp{Q, E) /dE is the normalized differential fragmentation function, i.e. the 
number of antiprotons between E and E + dE created by a parton jet of type 
j and energy Q. The fragmentation functions have been evaluated with the 
high-energy physics event generator PYTHIA/JETSET |18| based on the string 
fragmentation model. 

Once the spectrum of emitted antiprotons is known for a single PBH of given 
mass, the source term used for propagation can be obtained through 

j^^E) = rf^iMM^dM f^y'^ , (15) 

dEdtdV^ ' Jo dEdV '^'dMdVi \a{tform)J < Pm > 

where d'^n/dMdVi is the mass spectrum modified by Hawking evaporation until 
today, a(to) and a(tform) are the scale factors of the Universe nowadays and at 
the formation time t/orm (which is a function of the PBH mass), is the local 
halo density and < pm > is the mean matter density in the present Universe. 
The dilution factor, for t/orm ^ t^q, applies to all universes of interest. The 
last term converts the mean density into the local density under the reasonable 
assumption that the clustering of PBHs follows the main dark matter component. 
The quantity d^n/ dMdVi can be obtained through the mass loss rate which reads 
dM/dt = —a{M)/M'^ (by simple integration of the Hawking spectrum multiplied 
by the energy of the emitted quantum) where a{M) accounts for the available 
degrees of freedom at a given mass. With the assumption a{M) ^ const it leads 
to: 

Hence the spectrum nowadays is essentially identical to the initial one above 
M* = Sato ~ 5 X lO^'^g and proportional to below. 
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4 Propagation and source distribution 



The propagation of the antiprotons produced by PBHs in the Galaxy has been 



studied in the two zone diffusion model described in [|19], pOf. In this model, the 



geometry of the Milky- Way is a cylindrical box whose radial extension is i? = 20 
kpc from the galactic center, with a disk whose thickness is 2h = 200 pc and a 
diffusion halo whose extension is still subject to large uncertainties. 
The five parameters used in this model are: Kq, 6 (describing the diffusion coeffi- 
cient K[E) = KqPR^), the halo half height L, the convective velocity Vc and the 
Alfven velocity Va- They have been varied within a given range determined by 
an exhaustive and systematic study of cosmic ray nuclei data and chosen at 
their mean value. The same parameters used to study the antiproton flux from 
a scale-free unnormalised power spectrum in are used again in this analysis. 
The antiproton spectrum is affected by energy losses when p interact with the 
galactic interstellar matter and by energy gains when reacceleration occurs. These 
energy changes are described by an intricate integro-differential equation El 



where a source term qf^{E) was added, leading to the so-called tertiary compo- 
nent which corresponds to inelastic but non-annihilating reactions of p on inter- 
stellar matter. Performing Bessel transforms, all the quantities can be expanded 
over the orthogonal set of Bessel functions of zeroth order and the solution of the 
equation for antiprotons can be explicitely obtained |]19[. Thanks to this sophis- 
ticated model, it is no longer necessary to use phenomenological parameters, as 
in the pioneering work of MacGibbon & Carr |^, to account for the effect of the 
Galactic magnetic field. The propagation up to the Earth is naturally computed 
on the basis of well controlled and highly constrained physical processes instead 
of being described by a macroscopic parameter Tieak used to enhance the local 
flux. 

The spatial distribution of PBHs (normalized to the local density) was as- 
sumed to follow a usual spherically symetric isothermal profile where the core 
radius Rc has been fixed to 3.5 kpc and the centrogalactic distance of the solar 
system Rq to 8 kpc. Uncertainties on Rc and the consequences of a possible flat- 
ness have been shown to be irrelevant in The dark halo extends far beyond 
the diffusion halo whereas its core is grossly embedded within L. The sources 
located inside the dark matter halo but outside the magnetic halo were shown to 
have a negligible contribution. 



5 Experimental data and inflationary models 

The astrophysical parameters decribing the propagation within the galaxy being 
determined, for each set of initial parameters {^{Mpeak), Mp^ak,^,!) defining the 
mass spectrum given in section 1, a p-spectrum is computed. Fig. |l| gives the 
experimental data together with theoretical spectra for [3{Mpeak) = 5 x 10~^^ and 
l3{Mpeak) = 10-2^ while Mp^ak = M„ e = 1 and 7 = 0.35. The first curve is 
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PBH ontiproton spectro 




Figure 1: Experimental data from BESS95 (filled circles), BESS98 (circles), 
CAPRICE (triangles) and AMS (squares) superimposed with PBH and secondary 
spectra for j3{Mpeak) = 5 x 10~^^ (lower curve) and l3{Mpeak) = 10~^^(upper 
curves). In both cases, Mp^ak = M^, e = 1 and 7 = 0.35. 

in agreement with data whereas the second one clearly contradicts experimen- 
tal results and excludes such a PBH density. It should be emphasized that the 
computed spectra are not only due to primary antiprotons coming from PBHs 
evaporation but also to secondary antiprotons resulting from the spallation of 
cosmic rays on the interstellar matter. The method used to accurately take into 



account such secondaries is described in |^ and relies on a very detailed treat- 
ment of proton-nuclei and nuclei-nuclei interactions near threshold thanks to a 
fully partonic Monte-Carlo program. The uncertainties associated with the theo- 
retical description of cosmic-rays diffusion in the Galaxy (coming from degeneracy 
of the model with respect to several parameters, from nuclear cross sections and 
from a lack of measurements of some astrophysical quantities) are described in 
20] & pri and are taken into account in this work. To derive a reliable upper 



limit, and to account for asymmetric error bars in data, we define a generalized 



as 

i {(^i + \Qi)) 

where a*^^ and a'^^P^ {a^^~ and cr*^^'^") are the theoretical and experimental pos- 
itive (negative) uncertainties, $*'^(Qj) and are the theoretical and experi- 
mental antiproton fluxes at energy Qi. Requiring this to remain small enough, 
a statistically significant upper limit is obtained. 
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Figure 2: Maximum allowed value j3{Mpeak) as a function of Mp^ak with 7 = 0.35 
and e = 0.5, 1, 2. The gravitational constraint is computed consistently assum- 
ing critical collapse from a bumpy mass variance at all scales. The antiproton 
constraint is significantly stronger than the gravitational constraint in the region 
< M, < lOOM,. 



The maximum allowed values of (3{Mpeak) obtained by this method are dis- 
played in Fig. as a function of Mp^ak for e = 0.5, 1, 2 with 7 = 0.35. As ex- 
pected, the most stringent limit is obtained when Mmax = {i.e. eMp^ak = M^). 
The curve is clearly assymetric because the mass spectrum is exponentially su- 
pressed at when Mp^ak < whereas it decreases as a power law when 
Mpeak > M*- This constraint is significantly stronger than the gravitational one, 
the requirement flpBH,o < ^m,o, displayed on the right hand side of the plot. In 
order to constrain inflationary models producing a bump in the mass variance, 
one has to compute the values Mpeak and P{Mpeak)- These will depend on the 
parameters of the inflationary model considered and can be usually traced back 
to the microscopic lagrangian. A numerical computation of I3{M) must be per- 
formed for each model using spectra normalized on large scales with the COBE 
(CMB) data for given cosmological background parameters, e.g. f2A,o = 1 — ^m,o 
|[T3[| . In particular, in a fiat universe with JIa^o = 0.7, the mass variance at the 
PBH formation time is reduced by about 15% compared to a flat universe with 

Our results differ from those obtained in p| in several ways. First, more 
experimental data are now available with much smaller errors as measurements 



from BESS98, CAPRICE and AMS 22 where added to the first resuhs from 



BESS93 Then, a much more refined propagation model is used. This is a 
key point as all the uncertainties on the astrophysical parameters used to describe 
the convective, diffusive and nuclear processes occuring in the Galaxy are care- 
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fully constrained and taken into account. The resulting antiproton flux can vary 
by more than one order of magnitude between extreme acceptable astrophysical 
models, making this study extremely important for the reliability of the results. 
Finally, the upper limit on (3 obtained in this work relies on PBH formation by 
near-critical collapse around the mass scale set by the bump in the mass vari- 
ance. Hence, in contrast with results obtained in |Q] , a constraint is obtained here, 
using eq.@, for different masses covering nearly three orders of magnitude. In 
addition, this allows us to obtain a constraint in the space of the inflationary free 
parameters for a given relevant inflationary model using the accurate expressions 
(i, m in (0). 

To illustrate how inflationary models can be constrained, we use here a so- 
called BSI model for which the quantities Mp^ak and P{Mpeak) can be found 
numerically using the analytical expression for its primordial power spectrum. 
The quantity F{k) is fixed by two inflationary parameters p and kg and exhibits 
a jump with large oscillations in the vicinity of kg, and the relative power between 
large and small scales is given by (an analytical expression for F{k) and relevant 
figures can be found in [^,|^). This feature derives from a jump in the first 
derivative of the inflaton potential at the scale kg so that one of the slow-roll 



conditions is broken and the resulting spectrum is quite universal p4|. Using the 



formalism of Section 1 one finds fcpeafe, which must be distinguished from kg, as 
well as P{Mpeak)- Numerical calculations give Mg = M(tkJ ~ 1.6 Mpeak- We are 
interested in spectra with p < 1, corresponding to more power on small scales. In 
Fig. 3, the constraint on the inflationary parameter p is displayed as a function 
of Mg. In other words each point in the plane Mpeak, P{Mpeak) is translated into 
the corresponding point kg, p. As p decreases, the bump in anitk) and P{M) 
increases. The constant spectral index n (already excluded by recent CMB data) 
which would pass successfully the antiproton constraint corresponds to n ^ 1.32, 
only slightly less than ?7,=1.33, the value satisfying the gravitational constraint 
at Mg ~ M^, 1^. Indeed, as mentioned in the Introduction, a small change in n 
gives a large variation in /5(M^,). 



6 Discussion 

Several improvements of our work can be expected in the forthcoming years. On 
the theoretical side, a better understanding of possible QCD halos appearing near 
the event horizon of PBHs should slightly alter the expected antiprotons fluxes. 
The very same computation should also be performed for gamma-rays, following. 



e.g. 0, and compared to the previously obtained limit on /3 in P] and HlO 
Although essentially independent, the results are expected to be close to the ones 
obtained here. 



On the experimental side, the AMS experiment should provide extremely 



accurate data of the antiproton flux on a very wide energy range. It should also 
allow to probe different solar modulation states, leading to a better discrimination 
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Figure 3: The minimal value of the inflationary parameter p is shown in function 
of Ms = Mnitk^) together with the gravitational constraint (straight hues). For 
given values (e, JIa^q), the region under the corresponding curve is excluded by 
observations. The three solid curves at the bottom (e = 2, 1, 0.5, from the left 
to the right) are the current constraints for I^a q = 0.7, the upper solid curve 
corresponds to ^mfi = 1 and e = 1. The two dashed curves, both for (1, 0.7) 
show the improvement expected if no antideuteron will be found (the lower, resp. 
upper curve refers to AMS, resp. GAPS). 
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between the signal and the background [^. Finally, it will be sensitive to low 
energy antideuterons which could substantially improve the current upper limit 
on the PBH density. According to , if no antideuteron is found in three years 
of data, the limit on (3{Mpeak) will be improved by a factor of 6. Furthermore, the 
GAPS project if actually operated in the future, would improve the bound 
by a factor of 40. 
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